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Abstract 

Consider the zero set of the random power series f(z) = ^ a n z n with i.i.d. complex 
Gaussian coefficients a n . We show that these zeros form a determinantal process: more 
precisely, their joint intensity can be written as a minor of the Bergman kernel. We 
show that the number of zeros of / in a disk of radius r about the origin has the 
same distribution as the sum of independent {0, l}-valued random variables X^, where 
P(X k = 1) = r 2k . Mor eover, the set of absolute values of the zeros of / has the 
same distribution as the set {ul. ^} where the Uf. are i.i.d. random variables uniform 
in [0,1]. The repulsion between zeros can be studied via a dynamic version where 
the coefficients perform Brownian motion; we show that this dynamics is conformally 
invariant. 



1 Introduction 

Consider the random power series 

f V ( Z ) = J2 a nZ n (1) 



oo 



n=0 



where {a n } are independent standard complex Gaussian random variables (with density 
e~ zz /ir.) The radius of convergence of the series is a.s. 1, and the set of zeros forms a point 
proc ess Zw i n the unit disk U. Zeros of Gaussian power series have been studied starting 



with 



Offordl (1965), since these series are limits of random Gaussian polynomials. In the last 
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decade, physicists have introduced a new perspecti ve, by interpret in g the zer o s of a Gaussian 
polynomial as a gas of interacting particles, see iHannavl (|1996T ). iLebceua (J1999J) and the 
references ther ein. Much of the recent inte rest i n Gaussian analytic functions was spu rred 
by the papers lEdelman and Kostlan and iBleher. Shiffman. and Zelditchl (J2000). A 

fundamental property of is the invariance of its distribution un der Mobius transformations 
that preserve the unit disk; see Section 2 for an explanation, and lSodin and Tsirelsonl ((2003) 
for references. 

Our main new discovery is that the zeros Zu form a determinantal process, and this 
yields an explicit formula for the distribution of the number of zeros in a disk. Furthermore, 
we show that the process Z^ admits a conformally invariant evolution which elucidates the 
repulsion between zeros. 

Given a random function / and points Zi, . . . ,z n , let p e (zi, . . . ,z n ) denote the probability 
that for all 1 < i < n, there is a zero of / in the disk of radius e centered at Z{. The joint 
intensity of the point process of zeros of /, also known as the n-point correlation function, 
is defined by the limit 

p e (zi, ...,z n ) 



p(z 1 ,...,z n ) 



lim 



(2) 



when it exists; see f|10|) for a related integral formula. 

Theorem 1. The joint intensity of zeros for the i.i.d. Gaussian power series (QJ) in the unit 
disk exists, and satisfies 

1 



p(zi, ...,z n ) = n n det 



[1-ZizA 2 ^ 



(3) 



Thus the zero set of the i.i.d. series fn(z) is a dete rminantal process in U, governed 



by the Bergman kernel K\j(z,w) = 7r _1 (l — zw)~ 2 ; see ISoshnikovl ([2000al ) for a survey of 
determinantal processes. In particular, (J3J) extends the known fact that p(zi, zi) < p(zi)p(z 2 ) 
for all Z\,Z2 G U, i.e., the zeros are negatively correlated. In fact, Z^ is the only process 
of zeros of a Gaussian analytic function which is negatively correlated and has a Mobius 
invariant law; see Section |21 

The determinant formula for the joint intensity allows us to determine the distribution 
of the number of zeros of /u in a disk, and identify the law of the moduli of the zeros. 



Theorem 2. (i) The number N r = \Z\j D B r (0)\ of zeros of fu in the disk of Euclidean 
radius r about 0, satisfies 



E(l + s 



*N r 



r 2k s) 



(4) 



k=l 



2 



for all real s. Thus N r has the same distribution as XlfcLi -^k where {X^} is a sequence 
of independent {0, l}-valued random variables with P(X^ = 1) = r 2h . 

(ii) Moreover, the set of moduli {\z\ : fu(z) = 0} has the same law as the set {£7^ }, 
where {Uk\ ar & i-i-d. random variables uniform in [0, 1]. 

From Theorem El we readily obtain the asymptotics of the hole probability P(N r = 0). 
Furthermore, the infinite product in (0J) occurs in one of Euler's partition identities, see (|36|). 
and this connection yields part (ii) of the next corollary. 

Corollary 3. (i) Let h = 47rr 2 /(l — r 2 ), the hyperbolic area of B r (0). As r f 1, we have 

-irh + o(h)\ /-it 2 + o(l) 



P(iV r = 0) = exp(Z^LtiW) = exp(: 



12(1 -r) 



(ii) The binomial moments of N r equal 

'_/V \ r k(k+l) 



E 



k J (1 -r 2 )(l - r 4 )---(l - r 2fe )' 



(iii) The ratio (N r — fi r )/a r converges in law to standard normal as r ] 1, where 



2 9 

fi r = EN r = , and a 2 = Var N r = 



1.1 General domains 



The covariance structure E^/u(z)/u(w) j = (1 — zw) 1 equals 2n times the Szego kernel 
Su(z, w) = (27r) _1 (l — zw)~ l in the unit disk. The Szego kernel So{z,w) and the Bergman 



kernel Kd(z,w) are defined, and positive de finite 



a smooth boundary. (See the next section or Bell (1992) for information on the Szego and 



or an y bounded planar domain D with 



Bergman kernels.) For such domains we can consider the Gaussian analytic function fo{z) 
with covariance structure 2hSd in D (an explicit formula for fu is given in ()12j)). Recall 
that a Gaussian analytic function in D is a random analytic function / such that for any 
choice of z±, . . . , z n in D, the random variables f(zi), . . . , f(z n ) have complex Gaussian joint 
distribution. 

Corollary 4. Let D be a simply connected bounded planar domain, with a C°° smooth 
boundary. The joint intensity of zeros for the Gaussian analytic function fr> is given by the 
determinant of the Bergman kernel 

p(zi, ...,z n ) = det[K D (zi, Zj)]ij. 

Note that for simply connected domai ns as in th e corollary, the Bergman and Szego 



kernels satisfy K D (z,w) = 4tiSd(z,w) 2 , see iBelll (J1992J), Theorem 23.1. 
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1.2 The one-parameter family of Mobius-invariant zero sets 

For p > 0, let Z U)P denote the zero set of 

Ai^) = E( / ) a ^ n ' (5) 



n=0 



where {a n } are i.i.d. stand ard complex Gaussians. In p articular, fin has the same d istrib u 



tion as fjj. As explained in lSodin and Tsirelsonl ([20031 ) (see also Bl eher and Ridzall ((2002)) 



for any p > 0, the distribution of Z^ tP is invariant under Mobius transformations that pre- 
serve U. Moreover, these are the only zero sets of Gaussian analytic functions with this 
invariance property. However, only p = 1 yields a determinantal zero process. 

Taking n — 1 in Theorem ^ one recovers the well-known formula (1 — \z\ 2 )~ 2 /tt for the 
i ntens ity of Zjj. More generally, the intensity of Zu tP in U is p/[vr(l — |-z| 2 ) 2 ], see 



Sodin 



( 20001 ) . It follows that the expected number of zeros in a Borel set A C U is p/(47r) times 
the hyperbolic area A(A) = J A (ijnpp • (Integration is with respect to planar Lebesgue 
measure.) This can also be inferred from Proposition |S] below. In Sectional we prove the 
following law of large numbers. 

Proposition 5. Let p > 0, and suppose that {Ah}h>o is an increasing family of Borel sets 
in V, parameterized by hyperbolic area h = A(Ah). Then the number N(h) = \Zy n Ah\ of 
zeros of /u p in A h satisfies 

v N(h) p 

hm — - — = — a.s. 

fi,->oo h Atx 

1.3 Reconstruction of |/u,p| from its zeros 

Theorem 6. (i) Let p > 0. Consider the random function /u jP , and order its zero set Zu }P 
in increasing absolute value, as {zk}^ = i- Then 

oo 

l/u, P (0)| = c p ne P/(2fc) l^l a.s. (6) 
k=i 

where c p = e( p-7-7 ^/ 2 ' p~ p l 2 and 7 = lim n ^^^ =1 1 — log nj is Euler's constant. 

(ii) More generally, given ( G U, let {CaJ^Li be Z VjP , ordered in increasing hyperbolic 
distance from (. Then 

\M P (0\ = c P (i - ici 2 )- /2 n eP/(2fc) rSf ■ (7) 

fc=l 1 - ^ k 



Thus the analytic function fv,p{z) is determined by its zero set, up to multiplication by a 
constant of modulus 1. 

This theorem is proved in Section |3 
1.4 Dynamics 

In order to understand the negative correlations for zeros of /u, we consider a dynamic version 
of the zero set Zjj. Denote by Zy(t) the zero set of the power series Yl^Lo a n(t)z n , where the 
coefficients a n (t) are independent stationary complex Ornstein-Uhlenbeck processes; in other 
words, a n (t) = e~ f / 2 W n (e*), where {W n (-)} n >o are independent complex Brownian motions. 

A direct calculation gives that for the process Z v , the intensity ratio p(zi, Z2)/[p(z 1 )p(z 2 )} 
is strictly less then 1 and decreases to with as the hyperbolic distance between Z\ and z 2 
tends to 0. This repulsion suggests that when two zeros get close, there is a drift in their 
motion that pushes them apart. However, this is not the case. Instead, we have the following. 

Theorem 7. Consider the process of zeros {Zu(t)} in the unit disk and condition on the 
event that at time t = there is a zero at the origin, i.e., G Zjj(0). The movement of this 
zero is then described by an SDE which at time t = has the form 

dz = adW, 

where W is complex Brownian motion, there is no drift term, and 

oo 

l/^ = l/u(0)| = Ci ne 1/(2fc) |^| a.s. 

k=2 

Heuristically, any zero of fu oscillates faster when there are other zeros nearby; this causes 
repulsion. 

Analogous processes Zj^it) can be defined in general domains, and we shall show in 
Section [7| that the family of processes Zu(t) is conformally invariant (no time change is 
needed). Theorem [7| can be extended to p ^ 1 as well. 

Conditioning to have a zero at a given location. It is important to note that the 
distribution of fy given that its value is zero at is different from the conditional distribution 
of fu given that its zero set has a point at 0. In particular, in the second case the conditional 
distribution of the coefficient a\ is not Gaussian. The reason for this is that the two ways of 
conditioning are defined by the limits as e — > of two different conditional distributions. In 
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the first case, we condition on |/u(0)| < e. In the second, we condition on fu having a zero 
in the disk B e (0) of radius e about 0; the latter conditioning affects the distribution of a\. 
See Lemma ITH1 in Section |3 



1.5 Hammer sley's formula 



The starting point of the proof of Theorem ^ is a general permanent formula for the joint 
intensity of zeros for Gaussian analyti c functions. A version for polynomials is due to 



Hammerslevl ([19561 ) and 



Friedman 



(1990). The pe rmanent form Q for Gaussian polynomials 



appears in the physics literature (|Hannavl (Jl996|)). Closely related formulae for correlations 
between zeros of random sections of a positive holomorphic line bundle over a compact 
complex manifold were established bv lBleher. Shiffman. and Zelditchl (J2000). 

The version we need is for Gaussian analytic functions, that are not necessarily polyno- 
mials. 

Proposition 8. Let f be a Gaussian analytic function in a planar domain D such that 
Ef(z) = for all z G D. Given points z\, . . . , z n G D, consider the matrices 

A = (e/(*)/(^J) ; B = (E/'(^)7(%)) ; C = (e/'^TW) • 
Assume that A is nonsingular. 

(i) The joint intensity for the zeros of f exists and satisfies 



p(z 1 ,...,z n ) 



E(\f(z 1 )..-f(z n )\ 2 f(zi 



f{z n 



Consequently, 



p(z u . . . , z n ) 



det(vrA) 
perm(C - BA^B*) 



det(vrA) 
;r when z 

denote the set of n-tuples of distinct zeros of f . Then for any Borel set B C D n we have 



(9) 



(ii) Assume that A = A(z\, . . . , z n ) is nonsingular when z\, . . . , z n G D are distinct. Let Z^ n 



E#(Bnz 



Zj^) dz\ . . . dz n 



(10) 



The proof of this proposition is given in Section |SJ 

In the derivation of Theorem ^ from proposition |H1 we use conformal invariance, the 
i.i.d. property of the coefficients, and the beautiful determinant-permanent identity ()26|) of 



Borchardt 



(185 
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Remarks on the liter ature . A nice introduction t o the theory of Gau s sian analyt i c func 



tions is given in 



jodinl 



(2000); for earlie r results, see Hammerslev (1956) 



Bogomolnv et all (|l992f ). iKostlanl Jl993). 



Edelman and Kostlan 
Shiffman and Zelditchl ( 2003f) 



Close to the topic of this paper are 
Determinantal processes are also being intensively studied, see 



( 199*5[ and 



iedman 



Ha nnav 



Sodin and Tsirelson 
|200oi 



199C 

199611 

2003) 



Soshniko\ 



Lcbccuf (199 



Theorem 
between zeros of 



[T] provides further evidence for the analogy, suggested by 
Gaussian polynomials and the Ginibre ensemble of eigenvalues of (non-hermitian) random 
matrices with i.i.d. Gaussian entries, which is known to be determinantal. 



2 Conformal invariance and preliminaries 

Complex Gaussian random variables. Recall that a standard complex Gaussian random 
variable a has density e~ zz /ir, expected value 0, and variance Eaa=l. A vector V of random 
variables has a complex Gaussian (joint) distribution if there is a determinstic vector Vq such 
that V — Vq is the image under a linear map of a vector of i.i.d. standard complex Gaussian 
random variables. 

If X, Y are real Gaussian random variables of mean zero, then X + iY is complex 
Gaussian if and only if X, Y are independent and have the same variance. A complex 
Gaussian random variable Z with EZ = satisfies EZ n = for any integer n > 1. 

The complex Gaussian power series. Recall the power series /u in jQ. A Borel- 
Cantelli argument shows that the radius of convergence of /u equals 1 a.s. Clearly, the joint 
distributions of fu(z k ) for any finite collection {z fc } are complex Gaussian, so the values of 
/u form a complex Gaussian ensemble. Since /u is continuous, its distribution is determined 
by the covariance structure 

00 

E(/u(*)Mio) = 5>^) n = ( x - 2S5 T 1 - ( n ) 

n=0 

The right hand side is 2ir times the Szego kernel in the unit disk; it suggests a natural way 
to generalize the power series /u. 

The Szego kernel. Let D be a bounded planar domain with a C°° smooth boundary (the 
regularity assumption can be weakened). Consider the set of complex analytic functions 
in D which extend continuously to the boundary 3D. The classical Hardy space H 2 (D) is 
given by the L 2 -closure of this set with respect to length measure on 3D. Every element of 
H 2 (D) can be identified with a unique analytic function in D via the Cauchy integral (see 



Belli Jl9_9_2J), Section 6). 
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Consider an orthonormal basis {ip n } n >o for H 2 (D); e.g. in the unit disk, take ip n (z) = A= 
for n > 0. Use i.i.d. complex Gaussians {a n }n>o to define the random analytic function 



f D (z) = V^TC^antpn^z) 



(12) 



n=0 



(cf. (6) in IShiffman and Zelditchl (j2003l n. The factor of y27r is included just to simplify 
formulas in the case where D is the unit disk. The covariance function of fo is given by 
2tiSd (z, w), where 

oo 

S D (z,w) = ^2?p n (z)il> n { 



w 



(13) 



ra=0 



is the Szego kernel in D. The Szego kernel Sd does not depend on the choice of orthonormal 
basis and is positive definite (i.e. for points Zj G D the matrix (So(zj, Zk))j± is positive 
definite). 

Let T : A — > D be a conformal homeomorphism between two bounded domains with C°° 
sm ooth bound ary. The derivative T' of the conformal map has a well-defined square root, 
Belli \mk p. 43. If {V>n}n>o is an orthonormal basis for H 2 (D), then {Vv ■ (ioT)}„> 



sec 



forms an orthonormal basis for H 2 (A). In particular, the Szego kernels satisfy 



S A (z,w) = T'{zfl 2 T{wyi 2 S D {T{z),T{w)). 



(14) 



When D is a simply connected domain, it follows from the transformation formula (fH^) that 
Sd does not vanish in the interior of D, so for arbitrary p > powers S D are defined. 

Let {r/ n } n >o be an orthonormal basis of the subspace of complex analytic functions in 
L 2 (D) with respect to Lebesgue area measure. The Bergman kernel 



K D (z,w) = ^rj n (z)rj n 



w 



n=0 



is independent of the basis chosen, see Neharil (1975), formula (132). 

The Szego random functions with parameter p. Recall the one-parameter family of 
Gaussian analytic functions /u. p defined in (jSJ). The binomial expansion yields that the 



covariance structure E( fu,p(z)fu,p(w) ) equals 



El 



-P 



n 



z w 



E 



[l-zw)- p =[2TtS v {z,w)y. (15) 

n=0 v 7 n=0 

The invariance of the distribution of Zu }P under Mobius transformations of the unit disk is 
a special case of the following. 
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Proposition 9. Let D be a bounded planar domain with a C°° boundary and let p > 0. 
Suppose that either (i) D is simply connected or (ii) p is an integer. Then there is a mean 
zero Gaussian analytic function f Dp in D with covariance structure 

e(/ D)P (z)7^H) = [2ttS d (z,w)} p for z,w G D. 

The zero set Zu, P of fz> iP has a conformally invariant distribution: if A is another bounded 
domain with a smooth boundary, and T : A — > D is a conformal homeomorphism, then 
T(Za iP ) has the same distribution as Zd, p - Moreover, the following two random functions 
have the same distribution: 

f A , p (z)=T'(zy/ 2 -(f D!P oT)(z). (16) 

We call the Gaussian analytic function fo, P described in the proposition the Szego ran- 
dom function with parameter p in D. 

Proof. Case (i): D is simply connected. Let \1/ : D — ► U be a conformal map onto U, and 
let {a n } be i.i.d. standard complex Gaussians. We claim that 

f{z) = ^\zyi^( n P ) aw (i7) 

n=0 ^ ' 

is a suitable candidate for fo, p - Indeed, repeating the calculation in (fTKj) . we find that 



E[f{z)f{w)) = [*'(z)*'H]^ 2 (l - V(z)Hf(w))-> 

= [V{z)W(w)y/ 2 ■ [2tts v (^(z), v(w))y . 

The last expression equals [27tSd{z, w)] p by the transformation formula ()14j) . Thus we may 
define fo, P by the right hand side of ([17)1 . If T : A — > D is a conformal homeomorphism, 
then $ o T is a conformal map from A to U, so (fTTj) and the chain rule give the equality in 
law Since T' does not have zeros in A, multiplying fu p o T by a power of T' does not 
change its zero set in A, and it follows that T(Z& tP ) and Zd, p have the same distribution. 

Case (ii): p is an integer. Let {i/j n }n>o be an orthonormal basis for H 2 (D). Use i.i.d. 
complex Gaussians {a ni> ... jn : ni, . . . ,n p > 0} to define the random analytic function 

f DtP {z) = (2tt)^ 2 Yl «m,-,nAi (z) • • • ^ np (z) ; (18) 

ni,...,n p >0 
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Sodinl (j2000n for convergence. A direct calculation shows that /d iP , thus defined, satisfies 



E(f D>p {z)fn, P H) = W Yl ^nAz)^nAw)--^n p (z)^ np (w) = [2nS D (z,w)Y 



ni,...,n p >0 



Figure 1: The translation invariant root process and a Poisson point process with the same 
intensity on the plane 



The transformation formula (JUjl implies that the two sides of (116)1 have the same covariance 
structure, \2tt S\(z, w)] p . This establishes (116(1 and completes the proof of the proposition. □ 

The general theory of Gaussian analytic functions implies that, up to multiplication by 
a deterministic analytic function, the random functions fu jP are the only Gaussian analytic 
functions wi th marginal intensity of zeros proportional to the hyperbolic area element. See 



Sodinl ([20001) for a proof. 



Similarly, the zeros of the Gaussian analytic function 

^) = E(^r) ^ 



n=0 

have distribution which is invariant under rotations and translations of the complex plane. 
Note that here p is a simple scale parameter: F CtP (z) = Fc t i{-Jpz). 

Letting p — > oo in the definition of /u >p , one recovers that the limit of the rescaled point 
processes p^ 2 Zn n is the zero s et of Fc i', this phenomenon and its generalizations have been 



studied bv lBleher et al.1 (j2000fl . 

Figure ^ shows a realization of the whole plane Gaussian zero process along with a 
Poisson point process of the same intensity. The orderliness of the zeros suggests that there 
is a local repulsion taking place. One gets similar pictures for the Szego random functions 
in the unit disk. The two-point intensity for zeros at the points r and is given by Q. The 
most revealing formula is the ratio p(0, r) / (p(O)p(r)), which shows how far the point process 
is from a Poisson point process, where this ratio is identically 1. For general p, with the 
notation s = 1 — r 2 this ratio equals 

1 + (p 2 - 2 p - 2) ( s p + s 2+2 p) + (p + l) 2 {s 2 p + s 2 +") - 2 p 2 (s 1+ p + s 1+2 p) + s 2+3 p 

O^f ( } 
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0.2 0.4 0.6 0.8 1 



0.2 0.4 0.6 0.8 1 



Figure 2: Relative intensity at z = and z = rasa function of r for p = 1, ~ ,|, j, § and for 
p = 1,4,9,16,25. 

and in the case p = 1 it simplifies to r 2 (2 — r 2 ). For every distance r, the correlation is 
minimal when p = 1 (see Figure |2j). For all values of p different from 1, for small distances 
zeros are negatively correlated, while for large distances the correlation is positive. 

When p = 1, the zeros are purely negatively correlated: this special phenomenon is 
explained by the determinantal form of the joint intensity. 
Remark. The Szego random function for p = 2 



n=0 



coincides with the n — » oo limit of the logarithmic derivative of t he characteristic function 
of a random n x n unitary matrix (see 



Diaconis and Evans (2001 



1). 



The analytic extension of white noise. Next, we show that up to the constant term, 
the power series fu has the same distribution as the analytic extension of white noise on the 
unit circle. Let B(-) be a standard real Brownian motion, and let 

r-2w 



l>ZTT 

u{z) = / Poi(z,e u )dB(t). 
Jo 



Here the integral with respect to B can be interpreted either as a stochastic integral, or as 
a Riemann-Stieltjes integral, using integration by parts and the smoothness of the Poisson 
kernel. Recall that the Poisson kernel 

Poi(^) = J-» (\^) = ^ (r^-= ~ 1 )= &tSv(*M - ^ 

27r \ 1 — zw J 2ir \ 1 — zw J 2tt 
has the kernel property 

"2tt 

Poif^, w) = I Poifz, e ft )Poi(e ft , w) dt . 



o 
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(This follows from the Poisson formula for harmonic functions, see lAhfforsl ()1978i h Section 
6.3). The white noise dB has the property that if f\, f 2 are smooth functions on an interval 
and f* = f fi(t) dB(t) then E/f/ 2 # = / f 1 (t)f 2 (t) dt. By this and the kernel property we get 
e(u(z)u(w)^ = Poi(z,w). Therefore if b is a standard real Gaussian independent of B(-), 
then 

u(z) = ^u{z) + b - (20) 

has covariance structure E[u(z)u(w)} = irdt.Su(z,w). Now if v, 1/ are mean complex 
Gaussians, then ${Evv' = 2E(9fo/9Rz/); thus (|TT|) implies that u has the same distribution as 

Remark. Similarly, since /u j2 is the derivative of Ylm=i a mZ m / \frn, the zero set Z U]2 can be 
interpreted as the set of saddle points of the random harmonic function 

00 

m=l 

in U. More generally, in any domain D, the zero set Zd.2 can be interpreted as the set of 
saddle points of the Gaussian free field (with free boundary conditions) restricted to harmonic 
functions. 

Joint moments of complex Gaussians. We will need the following known fact for the 
proofs of Hammersley's formula and Theorem ^ 

Fact 10. If Zi, . . . , Z n are mean jointly complex Gaussian random variables with covari- 
ance matrix A4jk = EZjZ k , then E(|Zi • • • Z n \ 2 ) = perm(A^). 

Proof. We will check that in general for jointly complex normal, mean random vari- 
ables Zj, Wj we have 

k 

E[Z t ■ ■ ■ Z n W, W n ] = J2U VZ 3 W a{j) = perm (?Z 3 W k ) jk , 

(7 j = l 



where the sum is over all permutations o 6 S n . (See the book of ISimonl ([19791 ) for a similar 
statement in the real case.) Both sides are linear in each Zj and Wj, and we may assume 
that the Zj, Wj are complex linear combinations of some finite i.i.d. standard complex 
Gaussian sequence {Vj}. The formula is proved by induction on the total number of nonzero 
coefficients that appear in the expression of the Zj and Wj in terms of the Vj. If the number 
of nonzero coefficients is more than one for one of Zj or Wj, then we may write that variable 
as a sum and use induction and linearity. If it is 1 or for all Zj, Wj, then the formula 
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is straightforward to verify; in fact, using independence it suffices to check that V = Vj 
has ~EV n V m = n\l{ m=n j. For n ^ m this follows from the fact that V has a rotationally 
symmetric distribution. Otherwise, \V\ 2n has the distribution of the nth power of a rate 1 
exponential random variable, so its expectation equals n\. □ 



3 A determinant formula in the i.i.d. case 

The goal of this section is to prove Theorem and Corollary The proof relies on the 
i.i.d. nature of the coefficients of / = fu, Mobius invariance, Hammersley's formula and 
Borchardt's identity (|2f))l. 
For (3 e U let 

T fi (z) = f==f (21) 
1 — pz 

denote a Mobius transformation fixing the unit disk, and define 

(1 - |/?| 2 ) 1/2 



l-Pz 



so that Ta(z) = T'Jz). 



Remark 11. Recall that for two jointly complex Gaussian random vectors X, Y, the distri- 
bution of Y given X = is the same as the distribution of Y with each entry projected to 
the ortho complement (in L 2 of the underlying probability space) of the subspace spanned 
by the components Xj of X. 

Proposition 12. Let f = fv and z\, . . . , z n G U. The distribution of the random function 

T zl (z)---T Zn (z)f(z) (22) 

is the same as the conditional distribution of f(z) given f{z\) = . . . = f(z n ) = 0. 

Proof. First consider n = 1. The assertion is clear for z\ = 0; here the i.i.d. property of 
the {ak} is crucial. More generally, for z\ = (3, by (fTBj) the random function / = tq • (/ oT^) 
has the same distribution as /. Since Tg(/3) = 0, from the formula 

oo 

f(z) = T /3 ( y z)^2a k (T (3 (z)) k 

it is clear that the distribution of Tp • f is identical to the conditional distribution of / 
given f{(3) = 0, whence the same must hold for / in place of /. The proposition for n > 1 
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follows by induction: to go from n to n + 1, we must condition (/ | f(z{) = . . . = f(z n ) = 0) 
on f(z n+ i) = 0. By the assumed identity for n points, this is equivalent to conditioning 
(T Z1 ■ ■ -T Zn ■ f)(z) on f(z n+ i) = 0. By Remark fTTl conditioning is a linear operator that 
commutes with multiplication by the deterministic functions T Zi . Applying the equality of 
distributions (f(z) \ f{z n+ i) = 0) = T Zn+1 (z) f (z) completes the proof. □ 

Fix Zi, . . . , z n G U and denote 

n 

T(z) = l[T Zj (z). (23) 
i=i 

Since T Zk (z k ) = and T Zk (z k ) = 1/(1 - z k z k ), we have 

n 1 

T(z k ) = T' Zk (z k ) ■ H T Zj (z k ) = J] T~~~= r ■ I I - **) (24) 



for each k < n. 



Corollary 13. Let f = fv and zi,...,z n G U. The conditional joint distribution of the 
random variables (f'(z k ) : k — 1, . . . ,n) given that f(zi) = . . . = f(z n ) = 0, is the same as 
the unconditional joint distribution of (T'(z k )f(z k ) : k = 1, . . . ,n) . 

Proof. The conditional distribution of / given that f{zj) = for 1 < j < n, is the 
same as the unconditional distribution of T ■ /. Since T(z k ) = 0, the derivative of T(z)f(z) 
at z = z k equals T'(z k )f(z k ). □ 

Consider the n x n matrices A and M, with entries 



A jk = Ef(zj)f(z k ) = (1 - ZjZ k 



-1 



M 7 - fc = (1 - z-Zk) 2 



By the classical Cauchy determinant formula, see iMuirl (J1923) p. 311, we have 



1 ZjZ k 



k,j ■> k<j 

Comparing this to (12411 . we see that 



det(A) = n|T'(^)l- (25) 



k=l 
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We will need the classical identity of iBorchardtl (J1855) (see also iMind (|1978h ): 

1 \ . / 1 \ . / 1 



perm 



■ det 

Xj + ykJ j)k \Xj + yk /j>k 



det 



(xj + y k ) 2 / hk 



(26) 



Setting Xj = z- 1 and y k = —z k and dividing both sides by FJ ■ Zj, gives that 

perm(A) det (A) = det(M). 



(27) 



We are finally ready to prove Theorem^ Corollary0]is a direct consequence of the conformal 
invariance in Proposition |U] and the way the Szego and Bergman kernels transform under 
conformal maps (see (JHJ)). 



Proof of Theorem [TJ Recall from (jSJ) that 



p(z 1 ,...,z n ) 



E(\f'(z 1 )---f'(z n )\*\f(z 1 ),...,f(z n ) = 0) 



7r n det(A) 

By Corollary [THJ the numerator of the right hand side of (|2"%|) equals 

E(|/(zi) • • • f(z n )\ 2 ) J] \V(z k )\ 2 = perm(A) det(A) 2 , 

k 

where the last equality uses the Gaussian moment formula of Fact [TUl and (|23jl. Thus 
p(z u ...,Zn) = 7r~ n perm(A) det (A) = n~ n det(M), 

by (J2Zj). 



(28) 



(29) 



□ 



4 The number of zeros of fu in a disk 

In this section we prove Theorem 121 and Corollary 01 In fact, the corollary only uses part (i) 
of the theorem, so we delay the proof of part (ii) to the end of the section. 



Lemma 14. Let r\ < . . . < r m , let Bj = B r .(0), and let Nj = #(Zu n Bj). Then 



EN 1 (N 2 - 1) • • • (N m - m + 1) = Ilt- 1 ) 1 ^ 1 ! 



(30) 



where the sum is over all permutations o of {1, . . . , m}, the product is over all cycles v of 
the permutation, \v\ is the length of v, and r u = Y\j eu rj. 
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Proof. Applying Proposition |H] (ii) to the set B 1 x • • • x B n we get 



EN 1 (N 2 -l)---(N m -m + l) = [ 



p(zi, . . . , z m ) dzi ■ ■ ■ dz n 

x---xB m 



= / det (K(zi, z,-)) dzi--- dz m . 

J BiX---xB m 

Expanding the determinant and exchanging sums and integrals we get a sum over all per- 
mutations of m elements: 

^sgn(cr) / K(zi,z ai )---K(z m ,z (Tm )dzi---dz m . 

fj JBxX-XBm 

For each permutation <r, the corresponding integral is a product over cycles v of o of 

I» = J K(zi, z 2 )K(z 2 , z 3 ) ■■ -K{z\ v \, zi) dzt •■ -dz\ v \. (31) 

where v is an ordered subset of {1, . . . , m} and each variable ranges over the disk B Vi of 
radius r Vi . The formula for the Bergman kernel gives 

oo 

K{z u z 2 ) = — — - = 7T- 1 Y,(n + l){z{z 2 ) n . 

Using this, we expand the product in (}3T|) into a sum of monomials in the variables {zj}™ =1 
and {z~j}™ =1 ; then we integrate term by term. Each monomial in which the exponents of Zj 
and ~Zj are different for some j integrates to 0. Thus in all remaining terms the exponents 
of Zj and Zj are the same. Since Zj always comes as part of a product ZjZj+±, the exponents 
of Zj and Zj + i have to be the same as well. This implies that in nonvanishing terms all 
exponents agree, and we are left with 

oo „ 

I u = 7r~ ^ / (n + l) 1 " 1 (z{zi ■ ■ ■ z\ u \z\ u \) n dzi ■ ■ ■ dz\ v \ . 

n=0 J 

Since 

(n+1) / \z\ 2n dz = 2n(n + l) [ s 2n+1 ds = nr 2n+2 , 

JB r (0) JO 

setting r v = Yl jev r s we get 

oo 



.2 



+1) = ' v 

1 - rl 



n=0 ~ 1 v 

Since sgn(cr) = 1)^' +1 , this completes the proof of the lemma. □ 
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Proof of Theorem [2] (i). Denote f3 k = E(^ r ), and let V be chosen uniformly at ran- 
dom from all permutations of {1, ... , k}. Let q = r 2 . Then by (|3U|). we have 



A=En(-i) M+i r4 



_ q\y\ ' 

where the product is over cycles y of P. Since the cycle containing 1 of V has length that 
is uniform on {1, ... , k}, and given that cycle, the other cycles form a uniform permutation 
on the rest of {1, . . . , k}, we get the recursion 

A^D-^r?^. ( 32 ) 

with /So = 1. Consider the generating function (3(s) = J2k>o PkS k - Multiplying by ks k 
and summing over k > 1, we get 

S /3'( S )=/3( S )( S V(s)), (33) 

where 

OO £ oo oo 

€=i y fe=i e=i 

We write in the form 

(log /?(*))' = #0, 

which we integrate to get 

oo oo / k \£ 00 

iog/?oo = - E E ^r 1 - = £ + ^ > 

fc=l £=1 ' fc=l 

where the constant term is zero as /?o = 0(0) — 1. Thus 

oo 

(3(s) = H(l + q k s). (34) 
fc=i 

Taking expectations of the identity 

fc=0 ^ ' 

gives 

oo 

E ( a + l)"- = E^ = /3( S ) , 

fc=0 

and this concludes the proof of (i). □ 
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Proof of Corollary El (i) Theorem H implies t hat P(N r = 0) = UkLii 1 ~ r2k ) and the 
asymptotics for the right hand side are classical, see iNewmanl ( 19981 ). p. 19. For the reader's 
convenience we indicate the argument. Let L = logP(A^ r = 0) = YlT=i l°g(l — r2k ) which 
we compare to the integral 

1 



log(l 



„2k\ 



dk 



—2 log r 



log(l 



dx. 



(35) 



2 log r 



Eoo 
n=l 



We have / + log(l — r 2 ) < L < I, so L = I + o(h). Since — log(l — e x ) 
the integral in ("3*""j) converges to — 7r 2 /6. But 21 ~* r = = ^ + o(h), and we get 



7r 2 /12+o(l) 
l-r 



ttt + o(h) , as claimed. 



(ii) One of Euler's partition identities (see Pak (2003), section 2.3.4) gives 



fk+i\ 



-q)---(l-q k ) 



(36) 



k=l k=0 

so the claim follows from (J34j) . 

(iii) The formulas for the mean and variance of N r follo w from the binomial moment for- 
mula (13*01 . Using the general central limit theorem due to ICostin and Lebowita ()1995f ) and 
Soshnikovl ()2000b| ) (p. 497) for determinantal processes, we get that as r — > 1, the normalized 



distribution of N r converges to standard normal, as required. Alternatively, the last claim 
can be easily verified by computing the asymptotics of the moment generating function di- 
rectly. Yet another way is to apply Lindeberg's triangular array central limit theorem to 
the representation of N r as the sum of independent random variables, as given in Theorem 
Hi). □ 



The joint distribution of the moduli of Z\ 



u 



Proof of part (ii) of Theorem [21 

The zero set of fn is determinantal with the Bergman kernel K(z,w). Let 



^ n— 1 

if n (^) = -V(i + i)( 

7T ' 



ZW) 



3=0 



Since K„.(z,w) — » K 



z, w) a s n — > oo uniformly on compact sets of U 2 , Proposition 3.10 of 



Shirai and Takahashil (|2003l ) yields that the determinantal point processes with kernels K n 
converge weakly, as n — * oo, to Z-g. Thus it suffices to prove that the set of absolute values 
{|Cj|}™=i of the n random points of the determinantal process with kernel K n , has the same 
law as {Uj 3 }^ =l where Uj are i.i.d. uniform on [0,1]. 
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For any Z\ , . . . , z n 



( 1 



7T 



Z2 



,Zl, Z\ ) KJz^Zn) \ 



Zri) Z\ Zfit Zn) ) 

z 'r x \ ( i o ... o \ / i 



\ 1 z n ... zl' 1 / 



V 



o o 



n 



J 



i \ 



vr 1 



3T 1 / 



Setting Zj = rje j we find that the joint intensity of {|Cj|}? = i, evaluated at {r 3 }^ =1 , equals 



(37) 



n 



7T 



det I lf„ , z fc ) ) n d0i • • • r n d0 n 

: / (e ^wnj^r 1 ) (e s g n (^nj =1 ^ _l ) *•! ^ • • • r n m« 



When we expand the sums, for u/t the integrand contains a factor of the form z?z 9 - with 



3 3 

Thus (jSZD 



p q, and therefore the integral vanishes. When a = r, we get (27r) n JI -r^ 
equals 



Now f7? J has density 2jx 3 in [0, 1]. Hence, the joint intensity of {Ui , . . . , Un n } is precisely 
(I38|) . This proves the theorem. 



Remark. The proof above is modeled after an argument of iKostlanl (|l992j) for the distribu- 
tion of the eigenvalues of a random complex Gaussian matrix. It is simpler than our original 
proof that relied on random permutations. 



5 Law of large numbers 

The goal of this section is to prove Proposition the law of large numbers for the number 
of zeros of /u iP . We will use the following lemma in the proof. 

Lemma 15. Let p be a Borel measure on a metric space S, and assume that all balls 
of the same radius have the same measure. Let ip : [0, oo) — > [0, oo) be a non-increasing 
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function. Let A C S be a Borel set, and let B = Br(x) be a ball centered at x G S with 
fi(A) = fi(Bn(x)). Then for all y G S 



?/>(dist(y, z)) dfi(z) < / ip(dist(x, z)) dfx(z). 
a Jb 

Proof. It suffices to check this claim for indicator functions ip{x) = li x < r \. In this case, 
the inequality reduces to 

/i(A n B r (y)) < n{B R (x) n B r (x)), 

which is clearly true both for r < R and for r > R. □ 

Proof of Proposition We have 

EN(h) = [ p(z)dz = —h. 
Ja ^ 

Let Q(z,w) = p(z,w)/(p(z)p(w)). Then by formula (JT5j) we have 

Q(0,w) - 1 < (7(1- \w\ 2 ) p . 

we denote the right hand side by i/;(0,w) and extend if) to U 2 so that it only depends on 
hyperbolic distance. 




A JA 



E(N(h)(N(h) - 1)) - (EN(h)) 2 = / / {p{z, w) - p(z)p(w)) dwdz 

(Q(z, w) — 1) p(w) dwp(z) dz 




A J A 



< ij)(z 1 w)p(w)dwp(z)dz 



'A J A 

Let .Br(O) be a ball with hyperbolic area h = 4nR 2 / (1 — R 2 ). Note that p(w)dw is constant 
times the hyperbolic area element, so we may use Lemma ITB1 to bound the inner integral by 

if)(0,w) p(w) dw — c I (1 — r 2 ) p (l — r 2 ) 2 rdr 
b r (o) Jo 

with S = 1 — R 2 . Thus we get 

Var N(h) = EN(h) + E(N(h)(N(h) - 1)) - (EN(h)) 2 < ^ + ^ C s p ~ 2 ds. (39) 

47T 87T J s 
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For p > 1 this is integrable, so Var N(h) < 0(h). For p < 1 we can bound the right hand side 
of by 0(hS p ~ l ) = 0{h 2 ~ p ). Thus in both cases, as well as when p — 1 (see Corollary 
Efiii)), we have 

Var N(h) < c{EN(h)) 2 ~ p 
with (3 = p A 1 > 0. For r? > 1//3, we find that 

iV^) - EN(k^) 



EN(k r >) 

satisfies EF fc 2 = Oik^), whence E^2 k Y^ < oo, so Yfc — > a.s. Monotonicity and interpo- 
lation now give the desired result. □ 



6 Reconstructing the function from its zeros 

The goal of this section is to prove Theorem El The main step in the proof is the following. 
Proposition 16. Let c' = e p / 2 ~ 7 / 2 . We have 



l/u,p(0)|=c;iim(l-r 2 )-^ 2 J] 
We first need a simple lemma. 



\z\ a.s. 



\z\<r 



Lemma 17. If X, Y are jointly complex Gaussian with variance 1, then for some absolute 
constant c we have 



Cov loglXl , loglKI 



< c 



E(XY) 



(40) 



Proof. Write Y = aX+/3Z, where X, Z are i.i.d. standard complex Gaussian variables, 
a = E(XY) and \a\ 2 + \j3\ 2 = 1. It clearly suffices to consider \a\ < 1/2. Since 

log \Y\ = log \/3Z\ + log |1 + aX/(/3Z)\, 



the inequality (j4*0|) reduces to 



Cov^log \X\ , log 



1 + 



aX 



(3Z 



< c\a\ 



(41) 



We will use the estimate 



E|log|l + A/Z|| < ci|A| for AG 



(42) 
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which can be verified by considering the positive and negative parts of log 1 1 + X/Z\ as 
follows. The positive part is handled using the numerical inequality log 1 1 + w\ < \w\ and 
the integrability of \Z\ . For the negative part, when |A| > 1, the density of |1 + X/Z\ is 
uniformly bounded in the disk of radius 1/2, so it remains to consider the case |A| < 1. Then 
Elog_ |1 + X/Z\ can be controlled by partitioning into the events 

G k = { e - k < \1+X/Z\ < e 1 ^}. 

Since P(G fc ) = 0(|A 2 |e- 2fe ), we get 



E(l Gfc log_|l + A/Z|) =0(A;|A 2 |e- 2fc ). 



Summing over k establishes ()42j). 
By conditioning on X, (142)1 yields 



E log IX] - lot 



1 + 



aX 



pz 



< ci 



E X log \X\ = c 2 



This bounds the first term (expectation of the product) in the covariance on the left hand 
side of (jJTJ). The second term (product of expectations) can be bounded by the same 
argument. □ 

Proof of Proposition [THl Assume that f = / U(/0 has no zeros at or on the circle 
of radius r. Then Jensen's formula ( Ahlforsl (1978), Section 5.3.1) gives 

i r 27r \z\ 

log|/(0)l = — / log|/(re ia )|^+ lo Sy' 

71 J ° zeZ,\z\<r ' 

where Z = Z VtP . Let \ f{re ia )\ 2 = af,Y, where 

a 2 = Var/(re ia ) = [2nS v (r, r)] p = (1 - r 2 )~ p 
and Y is an exponential random variable with mean 1. We have 



Elog|/(re* a )| 



logcx 2 + ElogF — plog(l — r 5 



7 



2 2 
where the second equality follows from the integral formula for Euler's constant 



J — — I e x logxdx. 
'o 



Introduce the notation 



g r (a)=\og\f(e* a r)\ + 
22 



plog(l - r 2 ) + 7 



so that the distribution of g r (a) does not depend on r and a, and Eg r (a) = 0. Let 

»2- 



1 

2^ 



g r (a)da. 



We first prove that L r — > a.s. over a suitable deterministic sequence r n j 1. We compute: 



Var L r = E 



(2vr) 2 



2tt /.2tt 



g r (a)g r (P) d/3da . 



'o ./o 

Since the above is absolutely integrable, we can exchange integral and expected value to get 

VarL r = — — / / E(g r (a)g r (p))df3da = — / E(# r (a)# r (0))cia. 
(2viT Jo Jo 27r Jo 

where the second equality follows from rotational invariance. By Lemma IT71 we have 

|E(/(re fa )7(0) 1-r 2 

E(g r (a)g r (0)) < 



Var /(r) 

Let e = 1 — r 2 < 1/2. Then for a G [0, it] we can bound 



1 — r 2 e 



2 e ia 



1 1 -rV Q | > 



which gives 



-Var L r < ^ 



c2eP 



e |ot| < e 

2r 2 sin f > f e < a< vr/2 

1 vr/2 < a < 7T, 



I /-t/2 I 

<e L - p + -J ada + n/2<l c'|loge| p=l 



p<l 



J ^,2gia|p 



c'e 1 "'' p > 1. 



By Chebyshev's inequality and the Borel-Cantelli lemma, this shows that, as r — > 1 over the 
sequence r n = 1 — n~^ ly ^ 1 ' p ' +s ' , we have a.s. L rn — > and 

\z\ plog(l — r 2 ) + 7 



E lo sv 



log 1/(0)1, 



z£Z,\z\<r 



or, exponentiating: 



-7/2/ 



1 — r 



2 W 2 TT li 



riT 1 — i/c»i- 



(43) 



\z\<r 



Since the product is monotone decreasing and the ratio (1 — r 2 )/ (1 — r 2 +1 ) converges to 1, 
it follows that the limit is the same over every sequence r n — > 1 a.s. 
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Finally, by the law of large numbers (Proposition EJ), the number of zeros N r in the ball 
of Euclidean radius r satisfies 

AT, = ^(1+ = a.,, (44) 

whence 

r Nr = exp(iV r log r) = e ~ p/2+o(1) a.s. 

Multiplying this with (J43|) yields the claim. □ 



Proof of Theorem (i) By the law of large numbers for N r (see also 

V \ = 7 + log N r + o(l) = 7 + log p - log(l - r 2 ) + o(l) . (45) 

|2ftl<»* 

Multiplying by p/2 and exponentiating, we get that 

TJ gp/(2fc) = e 7P/2 p p/2 (1 _ r 2)-p/2 (1 + o(1)) (46) 
kfc|<r 

In conjunction with Proposition 1161 this yields (|B|). 
(ii) Let / = / U)P and 

T(z) : " ^ 



By (Jill), / has the same law as 

7=(T'r /2 -(/oT). (47) 
Now T'(C) = (1 - ICI 2 ) -1 - Therefore 

oo 

7(0 = (1 - IC| 2 )- p/2 /(0) = c P Y[e"/^\z k \ a.s., 

k=l 

where {z^} are the zeros of / in increasing modulus. If T(^) = z^ then are the zeros 
of / in increasing hyperbolic distance from (. We conclude that 

oo 

7(o = Cp (i - icr)- p/2 n eP/(2fc) i T (^)i a - s - 

k=l 

□ 

For our study of the dynamics of zeros in the next section, we will need a reconstruction 
formula for \f{j p (0)\ when we condition on that G Zjj jP . 
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Lemma 18. Denote by Q e the event that the power series /u )P defined in |3J) has a zero 
in B e (0). As e — > 0, the conditional distribution of the coefficients a l5 a 2 ,a3,... given Q £ , 
converges to a product law where a 1 is rotationally symmetric, \ai\ has density 2r 3 e~ r2 , and 
a 2 , a 3 , . . . are standard complex Gaussian. 

Proof. Let ao, ai be i.i.d. standard complex normal random variables, and p > 0. 
Consider first the limiting distribution, as e — > 0, of a\ given that the equation ao+ai^fpz = 
has a root Z in B e (0). The limiting distribution must be rotationally symmetric, so it suffices 
to compute its radial part. If S — \ao\ 2 and T = \ai\ 2 , set U = p\Z\ 2 = S/T. The joint 
density of (S,T) is e _s ~*, so the joint density of (U,T) is e~ ut ~H. Thus as e — > 0, the 
conditional density of T given {7 < pe 2 converges to the conditional density given {7 = 0, 
that is te~*. This means that the conditional distribution of a x is not normal, rather, its 
radial part has density 2r 3 e~ r2 . 

We can now prove the lemma. The conditional density of the coefficients ai, a 2 , ■ ■ ■ given 
Q e , with respect to their original product law, is given by the ratio P(fi e | di, Q>2, ■ ■ -)/P{^e)- 
By Lemma E01 the limit of this ratio is not affected if we replace /u iP by its linearization 
a o + a \\fpz. This yields the statement of the lemma. □ 



Kakutani's absolute continuity criterion (see IWilliama (|l99lt ). Theorem 14.17) applied to 
the coefficients gives the following 

Lemma 19. The distributions of the random functions fv, p (z) and (fv,p(z) — a o)/ z ore 
mutually absolutely continuous. 

Remark 20. By Lemma ITHl conditioning on G Z\j jP amounts to setting oq = and chang- 
ing the distribution of a± in an absolutely continuous manner. Thus, by Lemma [T9*l given 
G Z-g,f> the distribution of the random function g(z) = fv tP (z)/z is absolutely continuous 
with respect to the distribution of the unconditioned fu !P (z). Hence we may apply Theorem 
El to g(z) and get that given G Z%j jP , if we order the other zeros of /u, p in increasing absolute 
value as {zk}^^, then 



l/u,p(°)l = = c P He>'M\z k \ a.s. (48) 



k=i 



7 Dynamics of zeros 

In order to understand the point process of zeros of /% it is useful to think of it as a stationary 
distribution of a time-homogeneous Markov process. 
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Define the complex Ornstein-Uhlenbeck process 

a(t) := e-^W{e% W(t) := B ^ 

v2 

where Bi, B2 are independent standard Brownian motions, and W(t) is complex Brownian 
motion scaled so that EW / (1)W / (1) = 1. The process {a(t)} is then stationary Markov with 
the standard complex Gaussian as its stationary distribution. First we consider the process 

00 

<ft(z)ft(z; D) = J2 W n (t)4> n (z), t > 

n=0 

where W n are independent complex Brownian motions and {ip n {z)} n >o is an orthonormal 
basis for H 2 (D). With t = e T we get the time- homogeneous process, 

00 

f T (z) = e~ T/2 Lp e r(z) = 22a n (T)ip n (z). 

n=0 

Then the entire process <pt(z) (and so f T {z)) is conformally covariant in the sense that if 
T : A — ► D is a conformal homeomorphism, then the process 

{ y/r$jv>t(T(z)) } t>0 

has the same distribution as (ft(z; A), t > 0. For this, by continuity, it suffices to check that 
the covariances agree. Indeed, for s < t, 

E(p s (z)<ft{w) = E<f s (z)<f s (w) 

so the problem is reduced to checking the equality of covariances for a fixed time, which has 
already been done in Proposition El 

It follows automatically that the process {Z£>(t)} of zeros of (p t is conformally invariant. 
To check that it is a Markov process, recall from Section |21 that Zpit) determines (pt up to 
a multiplicative constant of modulus 1. It is easy to check that ipt modulo such a constant 
is a Markov process; it follows that Zr>{t) is a Markov process as well. 

Remark 21. This argument works in the case p = 1. By replacing the i.i.d. coefficients 
a n in (0) with OU processes, it is possible to define a dynamic version of the p 7^ 1 case 
in the unit disk. The same argument as above shows that these are Markov processes with 
distribution invariant under Mobius transformations of U. 

Finally, we give an SDE description of the motion of zeros. Condition on starting at time 
1 with a zero at the origin. This implies that Wq(1) = 0, and by the Markov property all 
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the Wi are complex Brownian motions started from some initial distribution at time 1. For 
t in a small time interval (1, 1 + e) and for z in the neighborhood of 0, we have 



<p t (z) = W (t) + Wx{t)z + W 2 (t)z 2 + 0{z 3 ). 

If W / i(l)W / 2(l) 7^ 0, then the movement the root z t of (fit where Z\ = is described by the 
movement of the solution of the equation Wo(t) + W\{t)zt + W%{t)z1 = O(zf). Solving the 
quadratic gives 

Expanding the square root we get 

Since Wo(t) is complex, W^t) is a martingale, so there is no drift term. The noise term then 
has coefficient —1/Wi, so the movement of the zero at is described by the SDE (at t=l) 
dzt = — Wi(t)~ 1 dWo(t) or, rescaling time for the time-homogeneous version, for any r with 
a (V) = we get 

dz T = —da (r). (49) 

ai(r) 

The absence of drift in (|49J) can be understood as follows: in the neighborhood we are 
interested in, this solution z t will be an analytic function of the {VK n }, and therefore has no 
drift. 

For other values of p the same argument gives 

dz T = — 1 . . da Q (r). 

Of course, it is more informative to describe this movement in terms of the relationship 
to other zeros, as opposed to the coefficient a%. For this, we consider the reconstruction 
formula in Remark |2*UI which gives 

oo 



Wi\ = \f^ P m = c P U eP/m \ z ^ a - s - 



k=l 



This means that when there are many other zeros close to a zero, the noise term in its 
movement grows and it oscillates wildly. This produces a repulsion effect for zeros that we 
have already observed in the point process description. The equation (j4T?j) does not give a 
full description of the process as the noise terms for different zeros are correlated. 
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8 Hammersley's formula for Gaussian analytic func- 
tions 



A version of the fol 
is from 



owing theorem was proved by iHammerslevi ([19561 ) . The present version 



Friedman! ([19901 ). Appendix B. We say that a point process has integral joint 



intensity p if formula (jlUj) holds for its counting function N. 

Theorem 22. Let f n = a n z n + . . . + oo be a random polynomial, so that (do, • • • , a n ) has 
an absolutely continuous distribution with respect to Lebesgue measure on C n+1 . Then the 
integral joint intensity of zeros exists and equals 



lim(7re 

e-»0 v 



1 \f'{z l )---f'{z k )\ 2 da ---da r . 

f(zi)t=B e {0),i=l,...k 



(50) 



We also need the following consequence of Cauchy's integral formula. 



Fact 23. Let D be a bounded domain, and let B C D be a closed disk. Then for every 
m > there exist constants c m so that for every f analytic on D and every z e B the m th 
derivative satisfies \f^ m \z)\ < c m (j D \f(w)\ 2 dw) 1 ^ 2 . 

Proof. Cauchy's integral formula gives a uniform bound on f^ m \z) for z G B in terms 
of the L 1 -norm of the function on any circle in D about B. Integration yields a bound in 
terms of the L 1 norm on an annulus, which is bounded above by the L 2 norm onD. □ 

Next, we note some consequences of the Taylor expansion for Gaussian analytic functions. 

Lemma 24. Let f be a Gaussian analytic function defined on a domain D, and let B C D 
be a closed disk about Zq. Consider the partial sums of the Taylor series expansion about z : 

n 

fn(z) = ^a k (z - z ) k . 

k=0 



Then for all m > the m th derivative satisfies 

supE|^ m ) _/M| 2 

B 



as n — > oo . 



(51) 



Consequently, for all mi,m 2 > the covariance function of the derivatives of orders mi 
and m 2 of f n converges uniformly on B 2 to the covariance function of the corresponding 
derivatives of f . 
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Proof. Note that finite a.s. limits of jointly Gaussian random variables are are jointly 
Gaussian with finite variance. This implies that the derivative of a Gaussian analytic function 
/ is a Gaussian analytic function. Moreover, the Taylor series of / has jointly Gaussian 
coefficients. Consider the L? space of functions on the set X = Q x D with the product of P 
and Lebesgue area measure. Assume without loss of generality that B is centered at 0, and 
let f n (z) = YTj=Q a j z ^ ■ Since /„ is a projection of / in the space X to the subspace spanned 
by /o, . . . , f n , it follows that f n —>f in L 2 (X). By Fact 1221 we have 

Esup l/f) | 2 ^0 

B 

and therefore 

supE|/M -/(™)| 2 -> 0. 

B 

which implies the claim for the covariance functions. □ 

Corollary 25. Let g n be polynomial of degree n with i.i.d. standard Gaussian coefficients 
independent of f n . Then f n + g n /n\ approximates f in the sense of (fo"Tj) and for each n has 
coefficients with continuous joint density. 

We first show a preliminary version of Proposition |S] (i) will then follow from the integral 
formula and a general lemma about point processes. 

Proposition 26. Using the notation of Proposition^ assume that A = A(zx,...,Zk) is 
nonsingular when z%, . . . ,z n G D are distinct. Denote by N(A) the number of zeros of f in 
A. Then for any n disjoint bounded Borel subsets Ax, . . . , A n of D, we have 

n „ 

EY[N(Ai)= p(z 1} ... } z n )d Zl ...dz n} (52) 

i= l </Aix---xA n 

where the integrations are with respect to Lebesgue area measure. 

Proof. Case 1: / is a polynomial whose coefficients have joint density. This is a conse- 
quence of Hammersley's formula, Theorem [221 
Case 2: D is the unit disk or the whole plane. 

A Fubini argument implies that there is a dense set TZd of rectangles in D such that for 
R G 7Zd, almost surely / does not vanish on dR. 

It clearly suffices to show the claim when the Aj are disjoint elements oUZd. 

Let {/m}m>i denote the approximation of / by polynomials in Corollary [2H 
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For A e TZd, the argument principle implies that the number Nm(A) of zeros of fu i n 
A, converges a.s. to the number N(A) of zeros of / in A. 

As M varies, the random variables YYi=i Nhii-h-i) are uniformly integrable, as they are 
uniformly bounded in L 2 by Lemma 1271 

The covariance functions of J'm and f' M converge uniformly on each Aj x Aj to those of 
/ and /', whence the permanent-determinant formula (on the right of (JHJ)) for Jm converges 
uniformly on Ai x . . . x A n to the permanent-determinant formula for /. 

Applying formula (fTUj) to fu and letting M — > oo we see that it converges to the desired 
formula for /. 

Case 3: D is simply connected: the claim follows from the Riemann mapping theorem and 
Case 2. 

Case 4: A general domain D. It suffices to prove (fTU|) when each Aj is a closed square in 
D. Then we can find a simply connected domain D Q C D that contains all the Aj, and we 
apply Case 3. □ 

For simple point processes, the following Lemma implies Proposition (|SJ) (ii). 

Lemma 27. Consider a simple point process (a random subset) Z in a domain D with 
counting function AT (A) = #(Z D A). Suppose that for any disjoint Borel subsets A 1; . . . , A& 
of D, we have 

E]jN(A l )= [ p( Zl ,...,z k )d Zl ...dz k . (53) 

i= l JA 1 x-xA k 

Let Z Ak C Z k denote the set of k-tuples of distinct points. Then for any Borel set 
B C D k , we have 

E #(£? n Z Ak ) = [ p(z u ...,z k )d Zl ... dz k . (54) 



J B 

Proof. Note that both sides of (J54j) define a Borel measure on subsets B e D k ; thus it 
suffices to show the equivalence for the case when B = B% x ■ ■ • x B k is a product set. 
Consider a finite Borel partition V of D and denote 

k 

M k{V)= Yl #Wix--xQ fc nBn^)= UNiQiHB,) 

Ql,...,Q k Qi,...,Q k i=l 

where the sum is over ordered A;-tuples (Qi, . . . ,Q k ) of distinct elements of V. Then the 
hypothesis (|53*j) implies that 



EM k (V)= £ f 



p(z 1 ,...,z k )dz 1 ...dz k , (55) 

x---xQ k nB 
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where we sum over the same fc-tuples as above. Consider a refining sequence of partitions Vj 
of D where the maximal diameter of the elements of Vj tends to as j ; — ► oo. By definition, 
M k (V) counts the number of fc-tuples (z ly . . . ,z k ) G B where z ly . . . z k are points of Z in 
distinct elements of V . We deduce that 

M k (Vj) -> \BH Z Ak \ 

monotonically as j — > oo. Taking expectations and letting j — > oo yields (J3j)) . □ 
We now proceed to analyze the behavior of Gaussian analytic functions near their zeros. 

Lemma 28. Let f be a Gaussian analytic function in a domain D, and assume that for 
every z G D a.s. z is not a double zero of f . Then a.s. f has no double zeros. 

The Gaussian assumption is needed: consider (z — 7) 2 with 7 a continuous random variable. 

Proof. We may assume that there exists z G D such that W = f(z ) — E(/(z )) is 
not identically zero (otherwise there is nothing to show). Let g(z) = E(f(z)W)/~E\W\ 2 and 
h(z) = f(z) — Wg(z). Then g is a deterministic analytic function with g(z ) = 1, and h(-) 
is independent of W. By assumption all the zeros of g are not double zeros of /. Any other 
double zero of / would also be a double zero of if) = W + h/g. If h/g is a random constant, 
then if> a.s. has no zeros. Otherwise, a.s. ip' = {h/g)' has at most countably many zeros 
and they are a.s. not zeros of if) since W 7^ —{hj g)(£,j) a.s. by independence. □ 

Lemma 29. Let f be a Gaussian analytic function (not necessarily mean 0) with radius of 
convergence tq, and let M r be its maximum modulus over the closed disk of radius r < tq. 
There exists c, 7 > so that for all x > we have 



P(M r > x) < ce 



-yx 2 



Proof. Borell's Gaussian iso perimetric inequality (see iPollardl ( 2002 ): the inequality 



was also shown independently bv iTsirelson et al.l (|1976T 1) implies that for any collection of 
mean Gaussian variables with maximal standard deviation a, the max M of the collection 
satisfies 

P(M > median(M) + 6(t) < P(iV > 6) , (56) 

where N is standard normal. Now the median of M r is finite because M r < 00 a.s. Since 
the distribution of f(z) is continuous as a function of z, the maximal standard deviation a 
in the disk B r (0) is bounded. The mean version of the lemma follows by applying (|56[) to 
the real and imaginary parts separately, and the general version follows easily. □ 
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Lemma 30. Let f(z) = a® + a\z + . . . be a Gaussian analytic function. Assume that Oq is 
nonconstant. Let A e denote the event that the number of zeros of f(z) in the disk B t about 
0, differs from the number of zeros of h(z) = Oq + a\Z in B t . 

(i) For all 5 > there is c > (depending continuously on the mean and covariance functions 
of f) so that for all e > we have 

P(A e ) < ce 3 - 25 . 

(ii) P(A t | ai, a 2 , . . .) < Ce 3 , where C may depend on (ax, a 2 , . . .) but is finite almost surely. 

Proof, (i) By Rouche's theorem, if \h\ > \ f — h\ on the circle dB e , then / and h have 
the same number of zeros in B e . By Lemma I2H1 applied to (/ — h)/z 2 , we have 

P(max |/0) - h(z)\ > e 2 ~ 6 ) < c exp(-^e~ 25 ) < Cl e 3 . (57) 

Let be the annulus dB ait + B e 2-s, and consider the following events: 

D = {\a \ <2e 1 - 5 }, 
E = {\ ai \ <e~ s }, 

F = {mm \h(z)\<e 2 - 5 } = {-a E 0}. 

Note that PE C < c 2 € 3 and that EOF C -Do- Given D , the distribution of a® is approximately 
uniform on B 2e i-s (i.e., its conditional density is 0(e 2S ~ 2 )). Since P(E) tends to one as e — > 0, 
this implies that 

P(F n E I D ) = P(-a E 0, E\D ) < c 3 gl^gll < c ^/ e ^ = c ^ 
and therefore 

P(F) < P(F n E\D )P(D ) + PE C < c 4 ec 5 e 2 - 2S + c 2 e 3 < c 6 e 3 " 25 . 

Together with (|57ji. this gives the desired result. Since all our bounds depend continuously 
on the covariance function of /, we may choose c in a continuous manner, too. 
(ii) The argument used to bound P(F) in (i) also yields that 

P( min \h(z)\ < 2\a 2 \e 2 {%}j>i) < c 7 e 3 . 

An application of Rouche's Theorem concludes the proof. □ 

The following lemma relates the integral joint intensity to the pointwise (strong) version. 
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Lemma 31. Consider a simple point process on a domain D. Let Zj G D, j = 1, . . . ,n. 
Assume that there exists disjoint neighborhoods Dj of Zj and a 5 > so that the integral 
version of the joint intensity satisfies 

p(zi, . . . , z n , z*) < C2\Zj — z # \~ 2+s on Dx x • - • x D n x Dj for all j. 

Let Nj je denote the number of points in the ball of radius e about Zj. As e — > 0, we have 

P(N 1>e = . . . = N n>£ = 1) < E(iV lj£ • ■ ■ N n , e ) = P(iV lj£ = ... = N n , e = l) + o(e 2n ). (58) 

Proof. For nonnegative integers Nj we have 

n n n 

< n N d - JJ l^- = 1) < JVt • • • N n J2(N k ~ 1). (59) 

j=l j=l fc=l 

The left inequality is clear. For the right one, if for some k we have > 1 then the kth 
term on the right alone gives an upper bound. We apply (J59j) to the Nj i€ with small e and 
take expectations. We apply Lemma 1271 to the set B e (zi) x ■ • • x B e (z n ) x B e (zk) to get 

E(iV v ■ • • N n , e (N kie - 1)) = / p(w u . . . , w n+1 ) dwi--- dw n+1 = o(e 2n ). 

J B c (z!)x---xB € (z n )xB € (z k ) 

□ 

Lemma 32. Consider a Gaussian analytic function f in a domain D with mean zero ev- 
erywhere. Let Zi,...,z n G D, and assume that for each j, the random variables f'(zj), 
f[zx),..., f(z n ) are linearly independent. Then there exists neighborhoods Di of the z\ so 
that for each 1 < j < n, and (u>i, . . . , w n , w*) G D\ x • • • x D„ x Dj the integral version of 
the joint intensity is defined and satisfies 

p(wi, . . . ,w n ,w*) < c\vjj - wA 2 . 

Proof. By continuity, there exists bounded neighborhoods Di of the Z{ so that 

(i) for all G IJ-Dj and u>, G Di the random variables f'(w*) and f(wi), 1 < i < n are 
linearly independent and the determinant of the covariance matrix is bounded away 
from 0, 

(ii) for all distinct points G [jDi and to; G Di the random variables f(w*) and f(wi), 
1 < i < n are linearly independent. 
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3)i 



Part (ii) follows by considering the Gaussian analytic function (f(w*) — f(wj))/(w* — w 
which has a removable singularity at Wj. Taylor expansion at w implies that for w, z G |J D 
the conditional distribution of f'(w) given f(w) = f(z) = is Gaussian with variance 
bounded above by c\\z — w\ 2 . Therefore, 

E(|/ / W ■ ■ ■ /'K)/'K)| | f{wi) = ... = /K) = 0) < c 2 e 4 , 
where e is the distance between w* and the set {wi, . . . , w n }. Furthermore, 

d 2 

— — det Cov(/(wi), . . . , f(w n ), f(w*))\ w = det Cov(/(w 1 ), . . . , / («;,,)). 

and since the right hand side is bounded away from 0, we get 

| det Cov(/(wi), . . . f{w n ), /(w*))| > c 3 e 2 . 

Now the permanent-determinant formula implies the claim of the lemma. □ 

Proof of Proposition [HI 
Step 1. We first verify the equivalence of (jHJ) and ©• First note that when / has Gaussian 
coefficients, then the values and coefficients of / are jointly Gaussian. In particular, the 
expression (joTIj) equals 

E(|/'(*i) ■ • • /' (z n )\ 2 1 f( Zl ) = ... = f(z n ) = 0) 0(0, . . . , 0) 
= E(|?/'W-?/'W| J )ff(0,.,0) l 
where g(0, . . . , 0) = n~ n det(A)^ 1 is the density of the Gaussian vector X = at 
0, and V is the projection to the orthocomplement of the subspace spanned by the entries 
of X. Setting Y = (f {zjj) n . =v note that the projection (J — V) olY onto the subspace 
spanned by the entries of X is given by BA~ 1 X. For a column vector Z of mean zero, recall 
that Cov(Z) = E(ZZ*). Now 

Cov(VY) = Cov(Y - BA^X) = Cov(F) - Cov^A^X) = C- BA~ 1 B*. (61) 

This proves the equivalence of (jHJ) and Q- 

Step 2. Part (ii). By Lemma the point process of zeros is simple. Part (ii) follows from 
Proposition OBI and Lemma ETH 

Step 3. Part (i). Let F e denote the event that / has a zero in each of B e (z-\), . . . , B e (z n ). 
Since the function p(zi, . . . , z n ) is continuous, we have 

P(F e ) < EN(B e ( Zl ))---N(B e (z n )) (62) 

p(zi, . . . , Zn) dzi ■ ■ ■ dz n 



B e (zi)x-xB e (zn) 

p(z 1 ,...,z n )7T n e 2n + o(e 2n ) 
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If, for some j, the derivative f'(zj) is not linearly independent of {f{zA : 1 < i < n}, then 
p(zi, . . . , £ n ) = and the claim follows from (|62j) . Otherwise, the claim follows from (|62|) . 
and Lemmas l3*Tl and 13*21 □ 
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